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In this work, a numerical approach has been proposed to estimate the five 
single-diode circuit model physical parameters of photovoltaic generators 
from their experimental current-voltage characteristics. Linear least square 
method has been used to solve the system of three linear equations to express 
the shunt resistance, the saturation current and the photocurrent as a function 
of the series resistance and the ideality factor. Two key points have been used 
to solve the system of two nonlinear equations to extract values of series 
resistance and ideality factor. The advantage of the proposed method with 
respect of existing numerical techniques is that use only two key points of the 
experimental characteristic and need only two initial guesses and does not use 
any approximation. To evaluate the proposed method, three PV generators data 
have been used to compare the experimental and the theoretical curves. The 
application of the proposed method provides a good agreement with the 
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1. INTRODUCTION 

The photovoltaic (PV) generator is the basic element of the photovoltaic system. This element is made 
of semiconductor material with a p-n junction which has the property of, when exposed to sunlight, producing 
direct current electricity proportional to the solar irradiation. The behavior of PV generator has a non-linear 
current-voltage (I-V) characteristic. Various representations have been employed to describe the current- 
voltage (I-V) relationship in solar cells. In practice, the single diode model is the commonly equivalent circuit 
used to describe the I-V relationship. However, regardless of this model, there are important PV parameters 
that must be accurate and extracted for the simulation, design, performance, evaluation and control of PV 
systems. The five main parameters that describe the behavior of PV generator models are the generated 
photocurrent (Ipn), the saturation current (Io), the series resistance (Rs), shunt resistance or parallel resistance 
(Rp) and the diode ideality factor (n). The electric current produced by a solar cell depends on the intensity of 
the incident light and its properties. 

Several methods have been suggested in the literature to extract the five parameters. We can divide 
these methods into three categories: analytical, numerical and evolutionary methods [1]-[22]. Maouhoub has 
proposed an analytical method to determine the five physical parameters [15]. The key points have been used 
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to calculate the series resistance and the ideality factor. The linear least square technique has been used to 
determine to parallel resistance, the photo current and the saturation current. Ghani ef al. [16] provided a 
numerical method to extract the values of the series and shunt resistors using the Newton-Raphson algorithm 
and the Lambert W- function. A singularity problem, which is divided by zero, can occur if the initial conditions 
of the parameters are incorrectly chosen. Moreover, these methods require many assumptions to simplify the 
problem of extracting the five parameters. Zhang et al. [17] proposes a simple and efficient numerical method 
for extracting all the parameters of a solar cell from a single current-voltage (I-V) curve under constant 
illumination by reducing the parameters to three. The disadvantage of this approach is that it allows several 
simplifications by eliminating terms of the analytical equation of current (I-V). Hejri et al. [18] proposed a 
numerical procedure using three key points from manufacturer’s data sheets and by solving a system of four 
nonlinear equations; this method needs four suitable initial guesses. Recently, Stornelli et al. [19] proposed a 
new simplified method for the iterative estimation of the five PV module parameters. This method uses some 
approximations to determine the values of Rs and n. Yadir et al. [20] proposed another numerical method based 
on key points and solving a system of four nonlinear equations. El Achouby et al [21] proposed numerical 
method to extract the five physical parameters, operating at standard test conditions (STC). This method is 
based on the variation of the ideality factor and solving a system of four nonlinear equations. The problem of 
this iterative method is that needs suitable initial values of four physical parameters. 

In the present work, we propose a numerical approach to extract the five physical parameters based 
on the single diode model. In a first step, we formulate three analytical expressions of the parallel resistance 
Rp, the photocurrent I,, and the saturation current Ip using the linear least square method. These three 
expressions depend on the diode ideality factor n and the series resistance R;. Then, we extract these two 
parameters numerically by solving a nonlinear system of two transcendent equations at two key points, namely, 
a short circuit point and maximum power point. The proposed method has the advantage over other methods 
of employing only two key points of the experimental I-V curve and only two initial guesses and does not use 
any approximation and any computing of the slope at open-circuit voltage or at short-circuit current. To check 
the precision of our method, we extract the five physical parameters of the commercial RTC PV cell, the 
PWP201 PV module and the KC200GT PV module and we compare the theoretical current-voltage curve with 
the experimental one. We organize our paper into four sections. After an introduction, we present in the second 
section the one-diode equivalent circuit model and the proposed five-parameter extraction methods. In the third 
section, we use three photovoltaic generators experimental data to validate our proposed method and to 
compare with other methods; we evaluate the accuracy of the presented method by calculating the absolute 
error and the root mean square error. Finally, we close our paper with a conclusion. 


2. PROPOSED METHOD 
2.1. Single diode model 

In our study, we use the equivalent single diode circuit model to describe the behavior of the 
photovoltaic generator. This circuit is illustrated in the Figure 1. 


ph R 


Figure 1. Electrical circuit for single diode PV model 


The output current of the PV generator, with five parameters, is a function of the output voltage: 


I = Ipn ~ Ip (exp () -1) - (1) 


n.NsV th Rp 
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The five parameters of this model are: the shunt resistance Rp, the saturation current Io, the photo- 
current Ipn, the ideality factor n and the series resistance Rs. 
- Ng is the number of cells in series in the PV module (Ns = | for one solar cell), 
- Vin = kgT/q is the thermal voltage, 
- q is the electronic charge, 
- kg is the Boltzman's constant 
- T is the temperature in Kelvin. 


This (1) is a nonlinear transcendent equation. We can write the current I = f (V) by introducing the 
Lambert (2). 


_ Jotlpn _ _ Gp ve nVth LoamnheAw ( I9-Rs —exp (ee )) (2) 


14Rs.Gp  1+Rs.Gp R nV tn.(1+Rs-Gp) nV tn(1+Rs-Gp) 


2.2. Proposed five parameters estimation method 
2.2.1.Analytical expressions of Io, Ipn and Rp 
According to our previous works [11, 15], we rewrite the (1) as (3). 


I; = Inn — Ip. Ai — Gp. B; (3) 
where 
VimtRslim = 
fi = ane ( nNsVth ) : (4) 
Bi =Vim t+ Rslim 


Vi, mis the measured voltage of the PV generator and N is the number of measured points. In order to 
formulate the analytical expression of Ipn, Io and Rp, we use the linear least squares method based on minimizing 
the follow objective (5). 


N 


S=) Ue- hm? 
i=1 


(5) 


I, mis the measured current of the PV generator and Ii is the theoretical current given in the (1). After 
application of this method, we formulate the follow linear system of three equations [15]: 


N N N 
lpnN Ip) Ai GP) Bi =) Ieexp 
i=1 i=1 i=1 


N N N N 

=lpn Y Ait ly > A? + Gp ) Bir =— > Iexp At 
i=1 i=1 i=1 i=1 
N N N N 

—Ipn > Bi tly) BiAr+Gp ) BZ =— > lvexp Bs 
i=1 i=1 i=1 i=1 


(6) 
Where Gp = 1/ Rp 


The solution of the previous system, provides three analytical expressions of Ipn (n, Rs), Io (n, Rs) and 
Rp (n, Rs) as a function of ideality factor n and series resistance Rs. 


2.2.2.Numerical extraction of n and Rs 

In order to extract the series resistance R, and the ideality factor n, we use two transcendent equations. 
By setting (V, I) = (Vimpp, Impp) using (1), we establish a first transcendent equation linking maximum power 
point to two physical model parameters: 
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Vinny + Rs. 1 Vinny + Rs. 1 
Impp = Ipn( Rs) — Ip (n, Rs) (exp ( mp “| 1) mp PP 


n. Ny. Vin R,(m, Rs) 
(7) 


By setting V = 0 using (1), we find a second transcendent equation linking short-circuit current Isc to 
two physical model parameters: 


Rslsc 
R,(m, Rs) 


mas RsIsc 
Isc = Ipn(1, Rs) — Ip(m, Rs) | exp NVR 1 
REESE OE 


(8) 


Equations (7) and (8) are two implicit and nonlinear equations with two unknown n and Rs. To 
determine these two values, we solve numerically this system of nonlinear equations via Newton-Raphson 
method by applying the nonlinear solver “fsolve” in MATLAB environment. To solve this system, suitable 
values of initial guesses are required. The main advantage of the presented method is that use only two initial 
guesses and needs only two key points of the experimental curve and does not use any approximation or slope 
computing. The flowchart in Figure 2 shows the steps involved in carrying out the proposed numerical method. 


- Input experimental I=f(V) curve, 
- Input Vmpp, Impp. Isc, Voc, T and Ns 


Initialize Rs, n 


Solve the system of two nonlinear equations (7) and (8) using 
‘fsolve’ in MATLAB 


Get new values of n and Rs 
Solve the system of three linear equations (6) 
Get values of Iph, I0, and Rp 


- Generate theoretical I=f(V) curve using equation (2) 
- Calculate RMSE value 


End 


Figure 2. Flowchart of the proposed method 


3. RESULTS AND DISCUSSION 

In order to validate the proposed method, it is examined for three experimental case studies: the 
Silicon solar cell RTC France at 33 °C (case study 1), Silicon Module PWP201 at 45 °C (case study 2), and 
the KC200GT multi-crystal PV module operating under Standard Test Conditions STC (case study 3). To test 
the accuracy of our method and to measure the degree of precision of our mathematical model defined above 
for the extraction of the five parameters, two statistical indicators are selected: the Absolute Error (AE) and the 
Root Mean Square Error (RMSE) estimator defined as (9) and (10). 
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N 2 
1 2 
RMSE = nD Si 
t=1 
(9) 
AE = |I; — lim| (10) 


The Table 1 shows the physical characteristics extracted from the experimental I-V curve of the three 
photovoltaic generators RTC France at 33 °C, PWP201 at 45 °C and KC200GT at STC. The experimental I-V 
curves of these two cases are extracted from [1]. For the KC200GT panel, the experimental I-V data are 
extracted from the datasheet [23] under standard STC test conditions (irradiation level 1000 W /m2, AM 
spectrum 1.5 and cell temperature 25 °C). 


Table 1. Basic parameters for RTC France cell, PWP201 and KC200GT modules 


Parameter RTC (33 °C) PWP201 (45 °C) KC200GT (STC) 

Vmpp (V) 0.4590 12.6490 27.11 

Impp (A) 0.6755 0.9120 7.59 

Voc (V) 0.5727 16.7785 33.07 

Isc (A) 0.7605 1.0317 8.21 

Ns 1 36 54 

Ki (mA/°C) : 2 3.18 

Kv (mV/°C) 2 < -123 


3.1. Case study 1: RTC solar cell 

Referring to the numerical method proposed previously, using the data listed in Table 1, we extract 
the values of the five parameters for the RTC France solar cell. The values, thus obtained in comparison with 
other proposed numerical approaches [8], [17], [20] are summarized in Table 2. As can be seen, the approach 
proposed provides a minimum value of (RMSE). The presented method has the advantage of using only two 
initial guesses, reducing the research space to the two unknown parameters and does not use any approximation 
and any computing of the slope at open-circuit voltage or at short-circuit current compared to other methods 
that use four or five initial guesses. Several of these methods determine the values of the slope at open-circuit 
voltage or at short-circuit current with many approximations. 


Table 2. Comparison of the different methods for the RTC France solar cell at 33 °C 


n Tofu A] Iph(A) Rs(Q) Rp(Q) RMSE 
Proposed 1.4877 0.3477 0.7611 0.0349 42.9764 1.4938E-3 
Yadir 1.4905 0.3514 0.7611 0.0355 45.0472 1.5555E-3 
Zhang 1.48 0.3267 0.7607 0.036 60.24 2.0135E-3 
Cannizzaro 1.4899 0.3604 0.7611 0.0355 44.8997 6.3813E-3 


Figure 3 shows the experimental characteristic I-V and the theoretical curve using the five estimated 
parameters, for RTC France solar cell. Theoretical curves are determined using the (2) and using the five 
extracted parameters. As it can be ascertained, the theoretical curves are in very agreement with the 
experimental measurement. Figure 4 shows the plot of the absolute error for RTC France solar cell compared 
with different methods. It is distinctly indicated that the proposed method in this work has a small error. The 
absolute error does not exceed 0.003 A. 
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Figure 3. Experimental and theoretical I-V Figure 4. Plot of absolute error vs. voltage for RTC 
Characteristic for RTC solar cell at 33 °C solar cell at 33 °C 


3.2. Case study 2: PWP PV module 

The method previously studied and validated in study case 1, it is applied to the PWP201 module. The 
extracted values are compared to the work of [2], [7], [17], and are abridged in Table 3. As can be seen, the 
proposed numerical method provides a very low RMSE value compared to other methods. 


Table 3. Comparison of the different methods for the PWP201 module at 45 °C 


n TolwA] Ton(A) RQ) R,(Q) RMSE 

Proposed 1.324858 2.725395 1.033699 1.209192 625.327 2.4096E-3 
Bai J 1.18355 0.578317 1.03273 1.33419 1121.67 16.495E-3 
Chan 1.3936556 5.04906 1.03254 1.12766 1123 9.2398E-3 
Zhang 1.27394 1.597 1.0332 1.313 602.3 6.5819E-3 


Figure 5 shows the experimental characteristic I-V and the theoretical curve using the five estimated 
parameters, for PWP201 module. As we can find out, the theoretical curve is in very good concordance with 
the experimental measurement extracted from the technical data sheet. Figure 6 shows the absolute error plot 
for the PWP201 module compared to different methods. It is very clearly indicated that the proposed method 
in this work has a very low error value compared to other methods. The absolute error does not exceed 0.005A. 


——= Prposed 
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0.025 | r= Chan DSH 

—— Zhang 


© Experimental data 


0.2 - 0.005 
m= Model 
0- : : : 4 0. f » 1 
0 5 10 15 0 5 10 15 20 
v(v) v(V) 
Figure 5. Experimental and theoretical I-V Figure 6. Plot of absolute error vs. voltage for 
characteristic for PWP201 module at 45 °C PWP201 solar module at 45 °C 
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3.3. Case study 3: KC200GT PV module 

The values achieved by this method for KC200GT module are summarized in Table 4. The extracted 
values are compared with those of other research authors [15], [18], [19]. In Figure 7, we display the absolute 
error curve for the KC200GT module in relation to the several comparison methods. The proposed method has 
a very low error value and the absolute error does not exceed 0.1 A. 


Table 4. Companson of the different methods for the KC200GT module in STC 


TO[uA] __Iph(A) Rs(Q) Rp(Q) RMSE 
Proposed 1089 2.56 8.225 0.229 128.96 0.04 
Maouhoub 1.0758 2.12 8.211 0.233 132.88 0.07 
Stornelli 11 3.27 8.196 0.2185 164.2 0.07 
Hejri 1.34 171 8.21 0.220 951.93 0.15 


To study the effect of temperature T and irradiance level G on the I-V characteristic, we suppose that 
the influence of the temperature and irradiation on n, Rs and Rp are negligible [15]. The other two parameters 
Ipn and Ip are given by (11) and (12). 


[pn (T, G) = — =. (lpn sre +i = Tsre)) (1) 
Ih (T,G) — Gy. Voc (T, G) 


eo (RTOS)? 


I)(T,G) = 


(12) 


0.6 


—— Proposed 
05 == Maouhoub 

= Stornelli 

—— Hejri 


30 35 


—] 
wa 


10 15 20 
v(v) 


nN 
on 


Figure 7. Plot of absolute error vs. voltage for KC200GT at STC 


Where Ipn, stc is the photo-current in STC conditions, Gstc is the global irradiation under STC 
conditions (Gstc = 1000 W/m7?), Tstc is the temperature under STC conditions (Tstc = 25 °C) and K;j is the 
temperature coefficient of Ipn. The temperature and irradiation dependence of Vo is given by (13). 


G 
Voe(T,G) = Voesre + Ky(T = Tore) + 2. No.Venln (=) 
STC 


(13) 


Kj is the temperature coefficient of Voc and Voc, stc is open-circuit voltage in STC conditions. Figure 
8 presents the experimental I-V characteristic and the theoretical curve using the five estimated parameters, for 
the KC200GT PV modules. As can be viewed, the theoretical curve of the above-mentioned PV module is in 
very good agreement with the experimental measurement extracted from the data sheet under the conditions of 
constant temperature T=25 °C and at different irradiation levels. In Figure 9, we plot the experimental I-V 
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characteristic and the theoretical curve, for the KC200GT PV module for fixed irradiation 1000W/m? and at 
different temperatures. It is clear that the theoretical values of the model are in accord with the experimental 
measurements. 


G=400w/m? 


O Datasheet 
___ Model 


5 10 15 20 25 30 35 
Vv(V) 


Figure 8. Experimental and theoretical I-V characteristic for KC200GT module under varying irradiance 


G=1000w/m? 


Figure 9. Experimental and theoretical I-V characteristic for KC200GT module under varying temperature 


4. CONCLUSION 

This paper proposed a numerical method to extract and evaluate the physical parameters of 
photovoltaic generators. The equivalent single-diode circuit with five parameters is applied for modeling the 
electrical behavior of the PV generators. This technique is based on the analytical extraction of the saturation 
current, the parallel resistance and the photocurrent as a function of the ideality factor and series resistance 
using the linear least squares method. Then, the estimation of the series resistance and the ideality factor has 
been found by solving numerically a system of two nonlinear equations. The suggested numerical method 
reduces the number of the initial values to two and extract the five parameters without approximation. The 
application of this method for the RTC France solar cell, PWP201 and the KC200GT panels indicated a low 
error and a good agreement with the experimental data. 
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